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In the original derivation of appendixes A and B, the load 
fac tor n vas in appropriately included i n the definition, of T . 
' TEe' ' foLlovihg 'cEanges should -1 !)!' hoEeST"" 


Page 4, change 
to 


t aerodynamic time [ (pVBy/nm)t] 

T aerodynamic time [ (pVSvr/ni)t] 


Pages 15 to 2k, change 141 to p. 

Page 20 equation B5 first line, change 
Page 20 equation B5 third line, delete 
Page 20 equation B5 fourth line, delete n. 
Page 20 sixth line, change b x ** <\/Npn to b x 




wV’Ejj 


Figure k, change np to p. 


The section "Design Charts" should he interpreted in accordance with 
the implications of the above changes. In particular, the distinctions 
made between np and p are no longer pertinent. 
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Previous NACA reports have indicated that it is possible to develop 
angles of sideslip which may cause critical vertical-tail loads in abrupt 
rudder-fixed rolls from accelerated flight, but the reliability of methods 
for predicting these sideslip angles has not "been demonstrated. In this 
report expressions for calculating the sideslip angles in these maneuvers 
are derived from theoretical considerations, and numerical solutions are 
obtained, for a wide enough range of variables to permit construction of 
design charts. Comparison of the maximum sideslip angles obtained from 
the design charts and from flight tests with those obtained using a 
greatly simplified expression indicates sufficiently close agreement to 
warrant use of the simplified expression for first approximations in pre- 
dicting sideslip angles and vertical— tail loads occurring in rolling pull- 
out maneuvers for conventional ailerons. An approximate method for treat- 
ing cases of nonlinear directional-stability characteristics is presented 
which gives reasonably good results. The vertical-tail loads measured on 
one airplane in rolling pull-out maneuvers corresponded closely with those 
calculated by the simplest methods when the actual sideslip angles attainod 
were applied. 


INTRODUCTION 


Recently attention has been directed to the rolling pull-out maneuver 
as a condition in which crtical loads might be developed on the vertical 
tail through the attainment of large sideslip angles (reference l) . Sub- 
sequent flight tests have verified the fact that the vertical-tail loads 
in rolling pull-out maneuvers may exceed design loads based on other 
maneuvers. To indicate the order of magnitude of these loads approximate 
expressions were presented in reference 1 fpr estimating the maximum side- 
slip angles and maximum vertical— tail loads developed in this maneuver; 
it was indicated in reference 1, however, that flight values might exceed 
the values computed by these approximate expressions. Comparison of the 
sideslip angles determined in flight with those computed using the 
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approximate expression of reference 1 verified that- the approximate ex- 
pression underestimates the sideslip angles developed, in most cases 
by a factor of the order of 2. This result indicated that the usefulness 
of the approximate expression of reference 1 is limited to the purpose 
of that report; that is, to demonstrate the importance of the rolling 
pull-out maneuver. 

To provide information hotter suited to design purposes a more com- 
plete analysis has been made of the rolling-pull out maneuver. In the 
analysis a simplified expression suitable for preliminary design is 
developed for predicting the sideslip angle resulting from the rolling 
pull-out maneuver. Design charts which may be utilized for more prociee 
computations are presented, and the effects of such factors as nonlinear 
directional-stability characteristics are discussed. Flight data aro 
presented and compared with the analytical results. 

The determination of vertical-tail loads in rolling pull-out maneuvers 
resolves itself essentially into the determination of the sideslip angles 
developed. This is demonstrated by the agreement shown in figure 1 be- 
tween vertical-hail loads determined in flight and those computed by the 
simplest methods using measured values of sideslip angle, with no regard 
for sidewash effects, differences in the dynamic pressure at the tail 
from free-stream dynamic pressure, or possible yawing velocities. For 
this reason the present report is devoted exclusively to the determination 
of the sideslip angles developed in rolling pull-outs. 


SYMBOLS 


The following symbols are used throughout this report; 

A aspect ratio (b w 2 /S v ) 

a real part of complex root 

b imaginary part of complex root 

b w wing span, feet 

g acceleration due to gravity, 32.2 feet per second per second 
I x moment of inertia of airplane about X-axis, slug-feet square 

Ig moment of inertia of airplane about Z-exis, slug -feet square 

I a l%/mb v 2 
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i 0 

radius of gyration at out X-axis, feet 
k z radius of gyration about Z-axis, feet 
Ly load on vertical tail, pounds 
tail length, feet 
m mass of airplane, slugs 

n normal acceleration divided by acceleration of gravity 
p rate of roll, radians per unit aerodynamic time 

p rate of roll, radians per second 

q free-stream dynamic pressure, pounds per square foot ^ 

q^ dynamio pressure at tail, pounds per square foot 

r rate of yaw, radians per unit aerodynamic time 

r rate of yaw, radians per second 

S w wing area, square feet 

S-fc vertical— ball area, square feet 

s operational parameter 

t time, seconds 

V velocity of airplane along flight path, feet per second 
v component of. flight velocity along Y-axis, feet per second 
¥ weight of airplane, pounds 

(3 angle of sideslip (positive when right wing is forward), radians 
(3° angle of sideslip, degrees 

5 r rudder deflection, degrees 

9 angle between horizontal, plane and relative wind, radians 
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** 

damping factor (used in e^ 1 ^) 


X 

wing taper* ratio / ^5 °^ lor ^.. 

V root chord/ 


fi 

relative density coefficient (m/pS w b w ) 

p 

air density, slugs per cubic foot 

T 

aerodynamic time [ (pVS v/ /fpi)t3 


T r 

relative rudder effectiveness 

[(aw^)(-_ 7S _;j 

9 

angle of bank, radians 


* 

angle of yaw, radians 


L 

moment about X-axis, foot-pounds 

ar 

moment about Z-axis, foot-pounds 

Nt 

normal force on vertical tail, pounds 

Y 

force along Y-axis, pounds 


%t 

vertical tail normal force coefficient (N^/ %S^.) 

sc Nt /aa t 

slope of ourve of vertical-tail normal-force coefficient 
against angle of attack, per degree 

Cl 

lift coefficient (nW/qS w ) 


C Y 

lateral force coefficient (Y/gS w ) 

c 2 

rolling-moment coefficient •(L/q_S w b w ) 

M z 

increment of rolling-moment coefficient due to lateral-control 
deflection 

C n 

yawing-moment coefficient (N/qS^b^.) 

ZsC n 

increment of yawing moment coefficient due to lateral-control 
deflection 

°Yp 

SCy/dP 


cz p 

50j/3p 
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Cnp 3Cn./Sp 

C l r dCi/d(rb w /2V) 


C nr dC n /d(rb w /2V) 

C lp bC l /Z( V h v /2V) 

c np ' aon/aCpWsv) 

Y v (q.S v /mV)C y!3 

Lp (qS v TD v /iii3cx 2 )Cj (3 


% (q_s w b v M z s )c nj3 

Np (q.S w b w / m k z 2 )(^ w /2V)C n5) 

N r ( qS^/ml^ 2 ) (b v /2V)C„ 

I’d (qS^/mk/Kb^JC, 

P 

Lr ( qS^v/rntj; 2 ) (b w / 2V ) C J r 


P* 




parameters used in computing 



of C n against P 


B° ' for nonlinear curves 

max 


THEORETICAL ANALYSIS 


For the purposes of the theoretical analysis the rolling pull-out 
maneuver is. considered to consist of an abrupt aileron deflection in 
accelerated flight, the ■ rudder being held fixed. The normal accelera- 
tion and the aileron deflection are considered constant throughout the 
maneuver, and the angle 9 between the horizontal plane and the rel- 
ative wind is considered small enough so that cosS can be set equal to 
unity. These assumptions are conservative in that they will result in 
computed sideslip angles larger than those that would be obtained in 
actual flight maneuvers where a finite time is required to reach maximum 
normal Acceleration or maximum aileron deflection or where the normal 
acceleration Is unsteady or where the angle 9 Is large. The effect of 
differences in 9 on the magnitude of the maximum computed angle of side- 
slip will be small, but the effect of unsteady normal acceleration may be 
larger, though still conservative. 
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In the analysis the parameter (££i/Cig)(C i/Cnp 0 ) is substituted 
for the parameter (pb/2T)(Cj-yc n ^o) used in reference 1. 

The equations and methods used in the theoretical analysis are 
given in detail in appendixes A, B, and C. Appendix A gives the equa- 
tions for which numerical solutions are obtained in order to develop 
design charts. In appendix B a simplified expression is obtained for 
calculating the maximum sideslip angle developed in rolling pull-cuts. 
Appendix C describes an approximation made for the gravity component of 
force on the airplane which permits its inclusion in the equations of 
motion as a linear factor. 


Simplified Expression 

The theoretical analysis presented in detail in appendix A and 
appendix B leads to the results plotted in figure 2 from which the 
following simplified expression for the maximum sideslip angle devel- 
oped in rolling pull-out maneuvers is deduced: 




max 


to/Cj )(Cx,/C n(3 o) 


1/4 


( 1 ) 


In the derivation of this expression the value of C n ^ was assumed 

as Cl/i 6. This value Is about the mean of the values of Cp/18 and 
Cl/ 14 which would be deduced for aspect ratios of 6 and 10 and a taper 


ratio of 0.5 from reference 2. The relative insensitivity of this value 
to changes in both aspect ratio and taper ratio within current design 
limits is noteworthy. The values of Cjp presented In reference 2 

are based on lifting-line theory; refinements to these values based on 
lifting-surface theory are shown in reference 4. 


In the development of equation (l), it was aleo assumed that the 
adverse yawing-moment coefficient of the ailerons was given by 


£C n = (££i/G ' t p ) (C L /l6) 

This is the theoretical value for a wing of aspect ratio 8 and taper 
ratio 0.5 having ailerons extending over the outer 50 percent of tho span, 
as obtained by combining data in references 2 and 3. These references 
may also be used to determine values of £C n for other wing-aileron 
configurations . 
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In reference 3 only the induced yawing moment due to the allerone 
is considered. For large aileron deflections or for unconventional 
ailerons the profile drag effect may also he important. An expanded 
form of equation (l) 'which, may he used to account for small differences 
in AC n from that assumed for equation (l) is given hy equation (2). 

/ci p n .i. x (2) 

(Ml/Cl^iQxJC^) Cl 8 


For reasons discussed in more detail later, the validity of equation (2) 
decreases as the value of dC n departs from (££j/ c ?p ) . 

The sum of the values of and £C n used in deriving equation 

(l) is equal to that used in deriving the equivalent expression given in 
reference 1. The value of the constant l/b given in equation (l) is, 
however, twice that obtained in reference 1, which indicates that the 
derivation of reference 1 which is "based on static conditions 1 b over- 
simplified. 

In the next section of this report. Design Charts, the results of 
a more exact analysis indicate that equations (l) and (2), while satis- 
factory for the preliminary design of airplanes with conventional arrange- 
ments, may "be greatly in error for airplanes with unconventional lateral- 
control devices such as spoilers. 


Design Charts 


In order to provide data suitable for design purposes, and to show 
by comparison the applicability of equations (l) and (2), a numerical 
analysis was made in which the maximum sideslip angle developed for 
each of several combinations of variables was determined. The equations 
of appendix A used for the analysis involve only minor assumptions and 
these are such as to result in Slightly larger compu'ted angles of side 1 -’ 
slip than would actually be obtained. 

The range of variables considered covers the limits of conventional 
design practice. The analysis was made for the conditions of the V— n 
diagram Bhown in figure 3 • Calculations were made for the curve of 
Cp, = 0*9 (curve A—B in fig. 3) and at a high— speed' point for n == 8 - 

(point C in fig. 3)* Results obtained from this analysis are considered 
equally applicable to the region within the boundary shown in figure 3 • 
Compressibility effects are not considered in the analysis. 
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Along the normal acceleration— velocity curve , values of nu of 
30, 75, and 120 were considered for a 0 L of 0.9 and of - 120 for a C L 

of 0.35. Hie value of p. for an airplane with a wing loading of 40 
pounds per square foot and a span of 40 feet at sea level is about 13. 
Variations in the other parameters such as, vertical— tail size, dihedral 
effect, moment of inertia about the airpl&ne Z— and Z— axes , and wing 
aspect ratio and taper ratio were considered either individually or in 
combination where it appeared advisable. The combinations of parameters 
used in these computations are given in table I. Since the analysis 
was carried out on a dimensionless basis, the velocity and normal accel- 
eration for any particular airplane configuration may be calculated from 
the expression 


V = 8.02 / n db y f ee t p e r second 

J cl 


n = 32.2 - P b w . (njj.) 

(W/8 V ) 

In cases where the oscillations were divergent the maximum value of the 
sideslip angle was considered to be that attained in the first peal. 

The results of the numerical analysis are presented in figures 4 
and 5 in a form that permits easy interpolation for design purposes. 

The curves of figure 4 cover the part of the V-n diagram which is 
limited by maximum lift coefficient (curve A-B of fig, 3)* In figure 
4(a) the variation of 0° m ax with (dCz/Cjp) (CL/C n pO) is presented 

for various values of C n pO, dC n , and nfi for a value of CipO » 
-0.0010; corresponding data for a value of Czp 0 = 0 are shown in 
figure 4(b) . 

Similar curves for very high speeds and high normal acceleration 
(point C of fig. 3) are shown in figure 5. For purposes of comparison, 
h max as calculated from equation (2) with £C n set equal to 
(dCj/Cj ) (G^/l6) is shown on all the curves. Also the results of 

ir 

applying equation (2) to the case of dC n =0 is indicated in figure 
4 for_ comparison with the corresponding curves obtained from the numerical 
analysis . 


The curves of figure 4 indicate that for preliminary estimates of 
sideslip angles and corresponding vortical-tail loads the use of equation 
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(2) for values of AC n around (ADi/Cjp) (C l/i 6), that is, equation (l) 

is satisfactory, the percentage error "being for most practical configu- 
rations of a relatively low magnitude and the direction of the error 
"being conservative except for arrangements having low dihedral effect 
and low directional stability . The deviations in the latter case are 
greatest for the lowest values of nu where, from the standpoint of 
vertical— tail loads, the importance of the deviations would be less, 
since low values of nu represent low values of normal acceleration and 
hence of Cj, which correspond to low values of 0. 

The agreement shown in figure 4 between the design charts and the 
curve representing equation ( 2 ) with AC n = 0 is poorer than the agree- 
ment shown with = (Cj/Cip) (0 l/i 6) . This poorer agreement results 

from the fact that equation (2) neglects a phase relationship that 
exists between the effects of and C n p. This phase relationship 

is properly accounted for only where ACu. * (Ci/C^p) (Cl/i 6) as in 
equation (l), so that equation (2) becomes less valid as it departs 
from equation (l). The varying discrepancies indicated in figure 4 
between the results of the numerical analysis and of the application 
of equations (l) and (2) may be used as an Indication of the dis- 
crepancies that will arise from the use in equation (2) of other values 
of AC n . 

Results of applying the numerical analysis to high values of nu 
and low values of Cj, which together correspond to high speeds and 
high accelerations are shown in figure 5 and indicate that for this 
condition the use of equation (l) is decidedly conservative for all 
configurations. This condition is not considered too important as re- 
gards vertical— tail loads because the maximum amount of aileron control 
is generally not applied at the highest speeds, with the result that the 
loads are not critical at the highest speeds. These curve b are included 
however, as an indication of the range of applicability of equation (l) . 

The effects of independent changes in several other variables that 
were consiiered in the analysis are indicated in figure 6. The results 
in figure 6(a) indicate that, for the changes in conf iguration assumed, 
the differences are of secondary order. Figure 6(b) shows that the 
rate of movement of the aileron control within the limits indicated has 
only a small effect on the maximum sideslip angles attained. 


Discussion of Nonlinear Characteristics 

The preceding analysis has been carried out assuming linear varia- 
tions of Cn. with 0° for all configurations. In practice, however, 
these curves as' well as those for other stability coefficients are 
frequently nonlinear. Accordingly, an analysis was made to develop 
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methods for handling nonlinear variations of C n with 0° that would 
permit use of the simplified equations (l) and (2) or the design charts 
of figure h. For this purpose numerical calculations were made of the 
maximum sideslip angles developed in rolling pull-outs, using the equa- 
tions of appendix A, hut modified by using' appropriate initial conditions, 
and for simplicity, hy using the angle of .hank q> instead of the approxi- 
mation of appendix C. For the calculations C n p° was assumed nonlinear, 
Cn r and Cy^ 0 were assumed to vary consistently with C^ 0 ,- and all 

other parameters of the airplane remained constant. The various curves of 
C p against 0° oovered hy the calculations are believed to encompass 
roughly the variations usually encountered in practice. The variations 
assumed are shown in figure 7 together with the results of the calculations 
presented as values of maximum sideslip angle attained for various applied 
rolling-moment coefficients AC! 2* The parameter A Ci was used instead of 

(ACj/C^p) CCy/Cnpo) in the ahscissa of figure 7 because for the nonlinear 
case no single value of C n po could logically he used in the latter term. 

The curves of figure 7 indicate that for the cases considered the 
variations of 0°max with AC j are consistent and may he predicted hy 
the following purely empirical method: 

1. Denote by (C n p°)x the elope of the curves of C n against 0° 

through 0=0, by (G n p 0 )a the slope of the curve of C n 
against 0° at values of 0° beyond the break in the curve 
and hy 0* the sideslip angle at which the break in the curve 
of C n against 0° occurs. 

2. Assuming each of the slopes (C n p°)x and (Cnp °)2 to exist 

separately through 0=0, compute the curves of 0°max 
against (ACz/Cjp) (Ci/Cnpo) from the design charts. 

3. Through 0 = 0 draw the curve of 0°max against ACj, corre- 

sponding to (Gnp°)i. Denote' this curve as line A. 

r“ -j 

k. Through 0° = 1.50* ) 1 - (C nf3 °)x/(Cnp °)2 • at AC 2 = 0, 

w 

draw the curve of 0°max against A Ci corresponding to 
(Cnp°) 2 . Denote this curve as line B. 

The final curve is composed then of line. A frem 0 = 0 to the in- 
tersection of lines A and B, and of line B - from the intersection on to 
higher values of 0°. The curves computed from this method for curves II, 
IV, and V are shown in figure 7 for comparison with those computed hy the 
numerical analysis. A reasonable fairing of the intersection of lines A 
and B may he applied for greater accuracy. 
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This analysis was made only for curves of C n against 3° which 
could he approximated by two straight lines. For cases in which this is 
not sufficient, or for cases in which extreme accuracy is desired, solu- 
tions may he obtained hy use of a differential analyzer or hy a step—' by- 
step integration as in reference 5 « ■ .. — 

The generality of the method presented and conclusions indicated hy 
the curves of figure 7 is not, of course, established by the few cases 
considered. The results do offer promise that with further analysis the 
conclusions will be verified or other rational simplifications will be 
developed. In assessing the value of the methods given here It is of 
interest to note that it gave good agreement with the maximum sideslip 
angles computed for the airplane of reference 5 by step-by— step methods. 


FLIGHT-TEST RESULTS 


Flight data which may be compared with the theoretical results pre- 
viously discussed have been obtained bn two airplane 0 , one of which was 
flown with two different vertical-tail configurations. YiewB of the air- 
planes tested are shown in figure 3. A typical time history of a roll 
out of a steady turn is given In figure 9. It will be noted in figure 9 
that the maximum value of the vertical-tail load occurs at the time of 
maximum sideslip. For airplane 1, Ames flight data obtained in aileron 
rolls were used, and for airplane 2 at configurations 1 and 2, Langley 
flight data on rolling pull-outs were used. For airplane 3, the maneu- 
vers were not made steadily enough to permit correlation with the design 
charts or with equation (l), the normal acceleration for most runs being 
less steady than the time history shown in figure 9» 


Comparison of Flight and Theoretical Data 

For the airplanes for which flight data were available, there were 
insufficient data to permit accurate estimation of Cjp or of C n ^ so 

that correlation could not Justifiably be made with the design charts 
presented in the preceding sections of this report. As an ind ication of 
the applicability of equation (l), however, the value of CnR was esti- 
mated by the method shown in table II. The resulting sideslip angles are 
compared with values obtained in flight tests in figure 10. As a matter 
of interest the values of sideslip angle computed from the approx ima te 
expression of reference 1, that is, 

go = CL (p-bw/gV) 


8 OCn/^°) 


( 3 ) 
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are also shown In figure 10. For simplicity the change in sideslip angle 
denoted by Aj3° is used in figure 10 instead of the absolute sideslip 
angle of p°. 

For airplane 1, excellent agreement is indicated between flight data 
and equation (l) and correspondingly poor agreement for equation (3) . 

(See fig. 10.) 

For airplane 2 with configuration 1 the comparison indicates reason- 
ably good agreement between flight values of P and values computed from 
equation (l). 

For airplane 2 with configuration 2, the agreement between flight 
data and equation (l) is less favorable. 

Although the data for airplane 3 were not steady enough to permit 
their inclusion in the correlation, it is of interest that when the 
maximum accelerations were used in the computations the values of side- 
slip angle were consistently larger than those obtained in flight. 

There are several factors entering into the foregoing comparison 
that would explain, at least partially, the disagreements noted and which 
should be considered in the Interpretation of all the comparisons. These 
factors, it will be noted, are essentially defects in the basic data and 
hence represent limitations in the application to these airplanes of the 
design chartB as well as equation (l). One of these factors is the value 
of C n p used in the approximate expression. The method used for deter- 
mining this value in the present case, noted in table II, Involves the 
estimation of the values of BCN-fc/dat and T r from a knowledge of 
geometric properties of the airplane and of the value of && r /dP as deter- 
mined from steady sideslips. The methods used for estimating the values 
of SCjf-t/dat and T r are based on wind-tunnel data ’(reference 6) and 
remain to be verified by flight tests. For airplanes that are already 
flying, a preferable method of determining C n g from flight tests is 
indicated in reference 7* 

In addition, the methods do not attempt to take into account ration- 
ally the possible nonlinearity of the curves of c n against p which 
are frequently found in practice. This factor is discussed at length in 
a preceding section of this report. In this connection It is significant 
that the curves of 5 r versus P in steady sideslip were less linear 
for configuration 2 than for conf iguration 1. of airplane 2, and the 
agreement between flight and computed values of sideslip angle was not 
so good for configuration 2 as for configuration 1. 
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A third source of error results from the use of the term 
(^Cj/Cip) (Cl/i 6) for the adverse yawing-moment coefficient of the ailer- 
ons. Aside from the small differences arising from differences in wing 
and aileron- configurations from that assumed, the theoretical analysis 
from which this value was obtained ■ (ref erence 2 combined with reference 3 ) 
accounts only for the induced drag and not for the profile drag due to 
aileron deflection Vhich may in some cases be of significant value. 


Vertical-Tail Loads 

For airplane 3 > the flight data were obtained at the Ames laboratory 
from simultaneous rudder-fixed pull-ups and rolls and from abrupt rudder- 
fixed rolxs from steady accelerated turns. Both maneuvers were basically 
a sudden application of ailerons in accelerated flight and no differentia- 
tion is made between the data for the two maneuvers. 

The maximum loads on the vertical tail as obtained from pressure- 
distribution measurements taken while performing these maneuvers are com- 
pared in figure 1 with those calculated using the expression 

Lv = <TtjSt 3° 00 

da t 


The values of 3° and q-j. used in the expression were flight values 
corresponding to the time at which the loads were obtained, and no allow- 
ance was made for the effects of sidewa3h as discussed in reference 8, 
and q.t was assumed equal to a. However, the data were corrected for 
the load changes resulting from small inadvertent movements of the rudder. 
At the time of maximum sideslip angle the tail loads computed in this 
manner gave good agreement with the measured loads; at other times in the 
runs as indicated in the time history of figure 9 , effects of yawing ve- 
locity, and so forth, would have to be inc'luded to obtain correlation. 

The scatter indicated in figure 1 is partly accounted for by the accuracy 
with which the loads are determined (error estimated to be 5 to 15 percent, 
depending on the absolute magnitude of the load). It appears, therefore, 
that equation (V) is adequate for estimating vertical-tail loads when the 
correct sideslip angles are applied. 


CONCLUSIONS 


From a theoretical analysis of the motions of an airplane in a rudder- 
fixed, rolling pull-out maneuver and from comparison of the results of the 
analysis with flight data the following results have been obtained: 
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1. From numerical solutions to the theoretical equations design 
charts were developed for predicting the sideslip angles in rolling pull- 
outs for a wide rang9 of variables. 

2. A simplified expression for computing the maximum sideslip angles 
in rolling pull-outs was derived. The maximum sideslip angles computed by 
this expression were sufficiently close to those obtained from flight tests 
and from the design charts to warrant use of the expression for prelimi- 
nary estimates of the maximum sideslip angles and hence the maximum 
vertical— tail loads . 

3. An approximate method was developed for treating cases of non- 
linear directional-stability characteristics. From a limited comparison 
with results obtained from a numerical analysis of the theoretical ex- 
pressions, the approximate method appeared to be generally applicable. 

b. The vortical-tail loads in rolling pull-out maneuvers corre- 
sponded closely with those calculated by the simplest methods when the 
aotual sideslip angles attained were applied. 


Ames Aeronautical laboratory, 

National Advisory Committee 
Moffett Field, Calif., 


for Aeronautics, 
August 19b6. 



NACA TN No. 1122 


15 


APPENDIX A 

EQUATIONS FOR NUMERICAL ANALYSIS 


9) 


The solution to the linearized lateral equations of motion (reference 


( -p ^ + (-r Els. ) + ( - 2 $ aci Ji. ft*** 

A ' r * ' la J 


<1t ia 


C>, » C*h^ jl£ + Ch^ C*>7^ £fl 

f -v EsA + (&-*?*£.) + ( -tip Effi A , 4 J n j±. 

' ic 7 V dT ia 7 ^ / ia 7 ■ ia 


Cy-Cy^f. = - KlZ*? ^ ~ 


^ i wt»k* ’ \ i 



,+ &) tfMrs 


(f£p 5 6 Oy ) - 0 
N it : ' 


(AD 


was obtained by operational methods using the Laplacian operator, such 
that (reference 10, p. 2) 


f(s) <s / f(x) e -82 - dx 
°o 


si(s) -f (o) = / ^ e“ sx dx 
dT 


r 

'I 


(A2) 


The reduced equations, therefore, can be -written, 


J-o 


~ C D 


ClaV- 


-r __L -2*0 = 2ff 


£CiH 


Bi c 
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5 (-5s) + r(8-2s;)-2*fff^ = 

io x *c 

* (" £ ) + f + 5 ( 0 ~ 5° r s) ’ 



0 


(A3) 


provided the Initial rates of roll and yaw and the initial angles of hank, 
yaw, and sideslip are all zero. The solution to these equations in terras 
of angles of sideslip is given, therefore, hy 








K'-XX^) 

(1) (0) 




Cl) (B -£*(!). . 




(Ah) 


where the symbol L —1 stands for the Inverse of the operation indicated 
hy equation (A2). The reduction of this expression is normally obtained 
hy factoring the denominator of equation (A4) and making use of the ox- 
press ion 
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In the present case, however, the denominator represents a quartic 
for which there is no practical general factorization, so that either a 
numerical solution or simplifying assumptions are required to obtain 
quantitative results in terms of the derivatives. 

The design cherts presented as figures h and 5 were obtained irc*a 
numerical solutions using values of the derivatives presented in table I. 


APPENDIX B 

DERIVATION OP APPROXIMATE EXPRESSION 

Neglecting the terms and C^/Iq in finding the roots to the 

quartic, an assumption which is best for high-speed unaccelerated flight, 
equation (Ah) is written 


*■> 


B = < 


(-£)(' 

61 






o 1 


la"" 

K Bi a 


(2* 

ic •' 

^ si c 

(1) 

(0) 




(0) ( s - (- 2# 

ic ic 




( 0 ) 




This reduces to 


_/ a - HIE. _ faf ^ off + f- - 

~ \ ia ia s ) sic I ic s \ 

P = 


(Bl) 


°n r N 
ic ) 


aj 


sio 


(Be) 


[s - (Ci^/ia) 3 (s - a - ib) (s - a + ib) 
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where 



(Cny/ic) + 


and where, by further neglecti n g 


t (Cyp/s) — (Cjap/io)] /k 


as compared to 

S&CCnplVio), * ~ VeiCCnpli/lc) 


The part of equation (B<?) multiplying ^(Z£ n ii/si c ) reduces to 
— [s — (Cip/i a )] since “-(Cz r CL/ia) ^7 b© neglected as compared to 

C jp/i a . This part can be rewritten 

f_JL L_J) 

“ i 0 lbs \ s-a—lb e-a+ib / 


which, according to equation (A5), has the inverse transform (reference 10) 


^ P T 

~2djt — — / e 3 ^ sin bx dx 


i 0 b 


(B3) 


The part of equation (Bl) which multiplies 2^(£C^p/i a s) can bo 
rewritten 


2# 


/ Cn- 


ela 




1 - 




~Ia a 1 

s - (Cz p /i a ) 


r 


+ ibl } 


r 




l2ib[aj^iCZp/i a ) + ib] J L 2ib [a - {CiJ ia) - ibji ^ 


. s — a — ib 


s — a + ib 


which, according to equation (A?), has the inverse transform 



NACA TN No. 1122 


19 


2* 


*a 


^np ■< •, ig, C 2p 

-~if + i CL ~i CL c^i^ 


^Ipn^a 

C7 p (a2+2)2)i c + [( Clp/i a ) ~ a] 3 + b £ 


- , np 


CxPnr 


~^CL + 


2i , 


'2a 


+ b 2 


[ (C2 p /ia) - a] 3 + b s 


e 82 coe bx 


t C(0l P /la) ~ a 3 

+ — e 33 - sin bx 

L(C2 p /i a ) - a] a + b a 


(ClPn r /snDic) [a 2 -b 2 - a(CZp/i a )] /(a 2 + b 2 ) ^ 

e sin bx dx 

[C2 p /la) - a] 3 + b s 


(B4) 


By assuming now that C n p * (Gl/i 6) and the adverse yawing-moment 
coefficient of the ailerons £C n is equal to (Cl/i6) (££i/Ci^) , and 

changing the notation and the variable in order to simplify the results, j 
the sum of equations (B3) and (B4) can be written 


c 


h 




•xv 
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r 




^ ° z p Cn r 




< J%ic 


Ejb 


lpl?0 LpK(3i c / KErv 

"is“ + ~~ ( 1 " T r~' 

+ : *- 9 


£ ■ 


a i 2 + b x 2 '(%> - ax ) 2 + b x £ 


U 


3 1 f 

7 + -ic 
16 2 L 


/ 2a x - Lp \ 

1 “ %• ( T ) I 

v a x 2 + bx 2/ -J 


(Lp — a x ) 2 + bx 2 


LpNp 0 aiX cos b x x 


r ipbi( j§ -|lo) (% - ^-koXrtyi 


( ax 2 - bi 2 - ailp ^ 

V a x 2 •+ b x a / 


(Lp - a x ) 2 + bi J 


1 6 


.*0. 


Q aiX Bln 'b ;L x \ ?fyix 


(B5) 


where 


ax - |(T t + N r ) 


hi » a/n^ 


Equation (B 5 ) may be written 


• # v 
= f(Ep, i 0 ) 

(^2C L /C2 p C n(3 o) . 


As an indication of the magnitude resulting from this analysis,, the 
following approximate values were chosen: 
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N r = -0.126 - 0.050 Np 
Y v = -0.177 ~ 0.012 N p 
i c » 0.143 

and equation (B5) was plotted on figure 2 for various values of Np and 
Ip. The curves show that for the assumption mentioned, equation (B5) 
may he written with little error as 

3°max 1 

(£C iCl/C ZpC n g° ) k 

APPENDIX C 

APPROXIMATION FOR SINq) IN THEORETICAL ANALYSIS 

The assumption made in solving equation (Al) that cp is equal to 
slop is equivalent to replacing the sine curve with a straight line 
having the ‘same slope as the initial slope of the sine curve, and "becomes 
Increasingly erroneous as 9 "becomes greater. A hotter approximation 
may he obtained hy finding the slope of a straight 'line which has the 
same integrated effect as the sine curve. This relationship may ho ex- 
pressed mathematically hy 



where 3c represents the desired straight line slope and T is the time 
of maximum sideslip. 1 

The angle of hank will certainly he greater than zero in the re- 
gion considered and may he replaced hy same average value 9 , so that 
equation (Cl) can he written as 
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In order to~solve equation (C2), an iteration process is used. That is, 
equation (Al) is solved with the original substitution of <p for stafp 
to "determine the variation of <p with t ami the valuo of Tj • Those 
values are used, in equation (C2) and k is determined. This value of 
k is then multiplied into the term ■fjCifp of equation (Al) and equation 
(Al) is again solved, this time for p. This second iteration usually 
is sufficiently accurate for the evaluation of Pretax; but if a check 
solution for <p and T x shows that it is not sufficiently accurate, 
the process may be repeated. 
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"^Far configurations 11 and 12, ^Yo, and Cri p vara ccmtined only in the same ccmbinationn as In 

configurations 1, 2, and 3. 
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TABLE II.- YALUES ASSUMED FOR AERODYNAMIC PARAMETERS 
IN DETERMINATION OF DIRECTIONAL STABILITY 
3C n /SB OF TEST AIRPLANES 


— — ~ ■ 

Parameters 

Airplane 

1 

(total, 

two 

tails) 

Airplane 

2 

Configu- 
ration 1 

Airplane 

2 

Configu- 
ration 2 

Total vertical tail area, 

S t , sq ft. . . . . 

91.0 

26.58 

23.72 

Rudder area (aft hinge line), 
sq ft . 

33.6 

8.65 

8.30 

Balance area, sq ft 

9.56 

1.97 

1.96 

Height (center— line stabi- 
lizer to tip along hinge 
line), ft 


6.51 

5.20 

Height along hinge line, ft. . 

8.43 

— 

— 

Effective aspect ratio of 
vertical tail . . 

1.56 

2.47 

1.77 

H, ft 

28.70 

18.59 

18.59 

dCif t /dat 

.038 

.049 

.041 

T r 

.655 

.585 

.615 

d5 r /d£° \ 

.620 

1.060 

.420 

dC n /d0° 2 

.000984 

.001575 

. .000490 

q.t/i 

1.0 

1.0 

1.0 


"^In steady sideslips from flight data. 
a ac n /^° is computed from the e sprees ion 


<&n (SCNt/^o-tJTr Cq.t/4) s t H (dS r /dB°) 
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Figure 1.- Comparison of vertical-tail loads computed with 
measured values of (3 with vertical-tail loads 
measured in rolling pull-out maneuvers in flight. Airplane 3. 








Figure 2.- Variation of tne term a 0 /(AOiCt /Cj Cnoo) e*-i m,-** 

JUCfcJl/ % - - *p - **p * - * & «_* * 

with Sp for various values of Lp for an 
airplane in high-speed unaccelerated flight* ic ^ 

0.143, Kx = -0.1364-0.0500 Sh, T t “ -0.117-0. Olid Nr. 
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(ACl/CiJ(Ci/C„ fi .) 

Figure 5.— Variation of the raa 
(AO i /G i p ) ( 0 i / Ongo ) 
and Cjpo. Configurations 11,13 
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Fisrure 7.- Effect of nonlinear variations of O n with P° 
6 on the maximum sideslip angle. 
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FIGURE «.- TWO -VIEW DRAWINGS OF THE AIRPLANES 

TESTED IN FLIGHT 
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Figure 9.- Typical time history of a roll out of a steady turn of 
airplane 3. 
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Figure 10.- Comparison of calculated values of Ago \ 
maneuvers in flight on several airplane 
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values -measured in rolling pull-out- 
figurations. 







